function Hv = H( a, b )
%H Возвращает вектор фокусировки для нашей антенной решетки
%   Входные аргументы - азимут и угол места

[kx, ky, kz] = sph2cart(b, a, 1);
k = [kx; ky; kz];
r{1} = [  1/8, 0, 0];
r{2} = [ -1/8, 0, 0];
Nant = length(r);

phi = nan(Nant,1);
for i = 1:Nant
    phi(i) = 2*pi * r{i} * k;
end

Hv = exp(1i*phi);
end

